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We consider a honeycomb network built of quantum wires, with each node of the network having 
a Y-j unction of three wires with a ring through which flux can be inserted. The junctions are the 
basic circuit elements for the network, and they are characterized by 3 x 3 conductance tensors. The 
low energy stable fixed point tensor conductances result from quantum effects, and are determined 
by the strength of the interactions in each wire and the magnetic flux through the ring. We consider 
the limit where there is decoherence in the wires between any two nodes, and study the array as a 
network of classical 3-lead circuit elements whose characteristic conductance tensors are determined 
by the quantum fixed point. We show that this network has some remarkable transport properties 
in a range of interaction parameters: it has a Hall resistance quantized at Rxy — h/e^ ^ although the 
longitudinal resistivity is non- vanishing. We show that these results are robust against disorder, in 
this case non-homogeneous interaction parameters g for the different wires in the network. 

PACS numbers: 



I. INTRODUCTION 

The transport properties of junctions of quantum wires 
are of interest both seen from basic and applied perspec- 
tives. From the basic physics aspect, quantum wires pro- 
vide experimentally realizable ways for studying inter- 
acting electrons in one-dimensional geometries, and in 
particular junctions where 3 or more wires meet can dis- 
play rather rich behaviors. Theoretically, the problem of 
quantum wire junctions is related to dissipative quantum 
mechanics in two or higher dimensions, and to bound- 
ary conformal field theory.l^ It also has a mathematical 
connection to certain aspects of open string theory in 
a background magnetic field. ^ ^ From a practical view- 
point, junctions of quantum wires should serve as im- 
portant building blocks for the integration of quantum 
circuits, as they are the natural element to split electric 
signals and serve as interconnects. 

Junctions of quan tum wires have been the subject of 
many recent studies, ^J^^^^^^ which have uncovered many 
interesting transport properties as function of interac- 
tion strength. Quantum wires with few transport chan- 
nels, at low energies, can be described as Tomonaga- 
Luttinger liquids, characterized by a Luttinger parameter 
g which encodes the electron-electron inter act ions .l^^^^^ 
The transport properties of a given junction depends on 
the Luttinger parameters for each wire. At low energies, 
the conductance properties of the junctions of n wires 
are encoded in an n x n conductance tensor or matrix 
Gij that relate the incoming currents to the applied volt- 
ages on the wires via li = Gij Vj . At low voltages 
and low temperatures, the tensor takes universal forms 
dictated by the nature of the infrared stable fixed points 
in the renormalization group (RG) sense. These fixed 
points have been categorized for the case of Y-j unctions 
(n = 3) of spinlesJ^ and spinfu^ electrons as function of 



the interaction parameter g when all the wires are iden- 
tical, and more recently in the case when the wires are 
not identical and have different values gi^ 
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FIG. 1: (a) Scheme of a grid showing the flow of the current 
and the boundary conditions. External currents are fixed, as 
well as the potential on the node on the upper right corner, 
(b) Building block of the grid: junction of three quantum 
wires with a magnetic flux threading the ring. The Vi,2,3 are 
the voltages applied on each wire, and the /i,2,3 the currents 
arriving at the junction from each of the three wires. 



In this paper we investigate the transport properties 
of networks constructed using Y-j unctions of quantum 
wires as building blocks. Fig. [l^ depicts an example of a 
network shaped in the form of a rectangle, and Fig. 
shows the individual Y-j unctions used in each node. We 
consider a simplified model where the 3x3 conductance 
tensor for each Y-j unction is taken to be that dictated by 
the low energy quantum RG fixed point, but the trans- 
port is treated classically between any two junctions. The 
treatment is sensible if the segment of the wires between 
two junctions is large compared to the characteristic de- 
phasing length in the system. But the length scales of 
the junction itself, for example the size of a ring as shown 
in Fig. [T]3, should be smaller than the dephasing length 
so that the junction is treated quantum mechanically. 
The case when the full system is treated quantum me- 



2 



chanically is extremely difficult to analyze, because it is 
an interacting problem. For instance, a lattice version 
of the problem would essentially be an example of a two- 
dimensional interacting lattice model with a fermion sign 
problem. 

We find rather remarkable results for the transport 
characteristics of the network of Y-j unctions, even when 
the role of quantum mechanics is just to select the RG 
stable fixed point conductances of the elementary build- 
ing blocks. When the conductance is controlled by the 
chiral fixed points x± we find that the whole network 
behaves as a Hall bar, with a Hall resistance that is quan- 
tized to Rxy = ±/i/e^, like in the integer quantum Hall 
effect, with the sign given by the particular chirality of 
the fixed points x+ oi" However, the longitudinal 

resistivity pxx 7^ 0, unlike in the case of the quantized 
Hall effect where pxx vanishes. The quantization of Rxy 
is a manifestation of the universal fixed point conduc- 
tances. The chiral fixed points are stable for a range of 
Luttinger parameters 1 < ^ < 3, and which of x+ oi" 
X- is selected depends on the fiux threading the ring in 
the Y-junction.^^ The flux breaks time-reversal symme- 
try, but it does not need to be quantized at any given 
value; because of interactions, the conductance of the Y- 
j unction flows to flxed point values for a range of fluxes. 

The quantization of Rxy = ±/i/e^ for the network as 
a whole is independent of the value of g in the wires, 
as long as they are in the range of stability of the chi- 
ral flxed points. Moreover, we show that the quanti- 
zation Rxy = ±/i/e^ is stable against disorder in the 
wire parameters. Speciflcally, we show that the quan- 
tization of Rxy remains even when the values of g for 
different wires are not uniform but disordered, i.e., they 
are randomly distributed around some average value g 
with some spread Sg. 

The paper is organized as follows. In Sec.|n|we briefly 
review the results for the conductance characteristics of 
single quantum Y-j unctions, which are the elementary 
building blocks for the honeycomb wire- networks. In 



Sec. Ill we present analytical results from which one can 



understand the origin of the quantization of Rxy when 
the conductance tensor of each of the Y-j unctions in the 
network is associated to a chiral flxed point. In Sec.[TV|we 
present numerical studies conflrming the analytical flnd- 
ings by analyzing grids with different values of the inter- 
action parameter g, different geometries and sizes, and 
extrapolate these results to the thermodynamic limit. 
These numerical calculations are of much value for the 
next step, taken in Sec. [V| where we discuss the robust- 
ness of the quantization of the Hall resistance in the case 
when the wires each have different Luttinger parameters 
distributed randomly. The Appendix contains a detailed 
description of the numerical method to solve our network 
of Y-j unctions. 



II. SINGLE Y-JUNCTION AS ELEMENTARY 
CIRCUIT ELEMENT 

Each of these Y-j unctions in the network consists of 
three wires that are connected to a ring which can be 
threaded by a magnetic flux, as shown in Fig. This 
flux breaks time-reversal symmetry, and the currents in 
the junction will depend on the potential at its extremes 
and the magnetic flux inside the junction. 

The current- volt age response of each Y-j unction is de- 
termined by its conductance tensor Gij. Within linear 
response theory, the total current Ij flowing into the junc- 
tion from wire j is related to the voltage Vk applied to 
wire k by 



(2.1) 



where j, /c = 1,2,3. Two sum rules apply to the con- 
ductance tensor because of conservation of current and 
because the currents are unchanged if the voltages are all 
shifted by a constant: 



- ^ Gjk 







(2.2) 



The Gij reach universal values at low temperatures and 
low bias voltages. These universal values are dictated by 
the RG stable flxed point that is reached for given values 
of the Luttinger parameters in the wires. Here we shall 
focus on the case where all the three wires have the same 
parameter g. In Sec.|V|we will consider the more general 
case of network of wires where the three wires for each 
Y-j unction have different ^'s. 

When the three wires have the same g, the flxed point 
conductance tensor has a Z3 symmetry and takes the 
fornPl 



Gjk — 



Gs 



{3Sjk - 1) 



Ga 



-^jk 



(2.3) 



where e^j = Sij-i — Sij-^i with i-\-3 = i and we separate 
the symmetric and anti-symmetric components of the 
tensor, whose magnitudes are encoded in the scalar con- 
ductances Gs and Ga- Ga vanishes when time-reversal 
symmetry is not broken, for instance in the absence of 
magnetic flux through the ring. 

The flxed point values of Gs and Ga depend on the 
strength of electron-electron interactions, encoded in the 
Luttinger parameter g. We will focus on the chiral flxed 
points x±7 which are stable in the range 1 < < 3.^^ 
In the chiral cases, the conductances are given by Gs = 
G^ = X3^^ ~ ^gG^. (As we will show in this 

paper, this form of Gs and Ga at the chiral RG flxed 
points leads to the quantization of the Hall resistance for 
the network.) Thus the chiral conductance tensors are: 



Gy 



[(3(5jfe - 1) ±gejk\ ■ 



(2.4) 
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We shall work in units where the quantum of conductance 
/h is set to 1. 

The Y-j unctions are then assembled into a network as 
shown in Fig. [l^. We consider a regular hexagonal grid 
of Y-j unctions with 2c external connections on both the 
top and bottom sides and r on both the right and left 
side. Parametrized in such a way and with wires of unit 
length, the dimensions of the grid as a function of r and 
c are 

Lx = 6c 

L^ = V3(2r + l). (2.5) 

In this grid we shall fix the current flow along the x- 
axis from left to right and we shall fix the currents fiowing 
into the top and the bottom to zero, as shown in Fig. 
[l^. Given the conductance tensors at every node of the 
network, we compute the potentials and the currents on 
the links of the grid. The resistances and resistivities of 
the networks are studied for different orientations and 
systems sizes, and for different values of ^. In appendix 
[A|we present details of the method used to numerically 
compute the response of the networks. 



which follows from rotational symmetry and the orienta- 
tion that we have chosen for the currents. 
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FIG. 2: Unit cell of the hexagonal network. Currents are 
assumed to be positive when directed along the arrows in the 
wires. Dotted lines denote the boundary of the unit cell. The 
rectangular region ABCD shown is used for computing the 
resistances and resistivities of the network. 



III. ANALYTICAL RESULTS 

We will measure the longitudinal and transverse re- 
sponses in the framework of the classical Hall problem 
by injecting a transverse current along the x-axis and 
imposing a zero current boundary condition along the 
two edges that are parallel to the x-axis. This approach 
suggests that we solve for the potential in the bulk as a 
function of the external current. In other words we need 
to invert the fundamental equation (2.1) for / and V for 
each wire in the bulk. 

While the full network problem is not tractable analyt- 
ically, we can still gain some insight from a combination 
of analytics and heuristics. In particular we will be able 
to prove quantization of the transverse resistivity ana- 
lytically, even with some forms of disorder. Similarly we 
will derive the general form of the longitudinal resistiv- 
ity. We will confirm these results numerically in later 
sections. Let us start with the unit cell of the hexagonal 
lattice. There are two vertices (nodes) in each cell and 
current is directed along the bonds (wires) as shown in 
Fig. [2] Looking at the right-hand node first, the poten- 
tials on the external wires, V2 and F3, and the potential 
on the internal wire Vi are defined only up to an additive 
constant. This means that Eq. (2.1) is not invertible. 



However, by setting Vi — 0, or equivalently shifting all 
potentials in the two nodes by a constant Vi - ^Vj — Vi, 
the gauge is fixed and we obtain, using Eq. (2.4), the 
following: 



V3 



^1 



1 



2 1T5 
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(3.1) 



The solution in the left node is similar but with the per- 
mutation (^2,1/3) ^ {V{.VJ,) and {h^z) ^ -(/a, ^3), 



Now consider the potential gradient in the x— and 
directions. It is straightforward to derive the change 
in potential per unit cell, A14 and AF^, directly from 
Eq. (3.1 ) as follows: 



= ^3 - ^2 = ^ [2/1 - /2 - /s (/2 - ^3)] 

^yV = 1^2 - ^3 = ^ + /2 - /3] . (3.2) 

It is instructive to consider a simple case. We will gen- 
eralize this result below, but for now consider a uniform 
current in the bulk in the x— direction (or "armchair" 
configuration to borrow nomenclature from graphene). 
Each horizontal wire in each unit cell has a current 
1\ = I. By symmetry the other wires split the current 
equally: I2 = h = I2 = I3 = —1/2. This configu- 
ration leads to a particularly simple potential gradient: 
A^V = 31 /2g and AyV = ±1/2. 

The result for the resistances and resistivities are ap- 
parent after we account for the geometric factors. Con- 
sider the rectangular region ABCD in Fig. [2j with sides 
dAB — >/3 and dAD =2. In the transverse direction 
the width of the rectangle is twice the distance between 
the midpoint of the wires (with currents I2 and /a), and 
the voltage drop Vab = 2 AyV. The Hall resistance 
(which coincides with the Hall resistivity pxy) is therefore 
Rxy = Vab /I = 2 AyV/ 1 = ±1. In other words the Hall 
resistance is independent of g and quantized to unity! 

Similarly, in the longitudinal direction the length of 
the rectangle is 4/3 the distance between the midpoint of 
the wires (with currents I2 and /s), and Vad =4/3 A^^ V. 
The longitudinal resistance is Rxx = ^/3 A^V / 1 — 2/g. 
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There is an additional geometric factor in the longitu- 
dinal resistivity given by p^x = {dAB /dAo) Rxx, and it 
is thus given by pxx = V^/g- Hence the resistivity is 
non-zero and there is dissipation unlike in the standard 
quantum Hall effect. 

Had we used an alternate ("zigzag") configuration 
where the transverse current is zero /i = and the uni- 
form current is in the direction, = = —I2 = 
—I2 = /, we would have found a similar result, i.e., 
that the resistance in the x— direction is quantized to 
Rxy = ±1 while the resistivity in the y— direction is 

Pyy V^/g- 

We find this result both unexpected and remarkable. 
By taking the classical conductivity limit for each wire we 
have allowed decoherence along the wires. However we 
have preserved the quantum coherence on each vertex, as 
the chiral relation Eq. (2.4) is by nature a consequence 
of quantum scattering. Nonetheless even after relaxing 
a portion of the coherence, some element of quantization 
in the thermodynamic limit has survived in the form of 
an integer quantized Hall resistivity. On the other hand, 
decoherence has destroyed the zero longitudinal resistiv- 
ity of the quantum Hall effect, and so we are left with 
a hybrid quantum-classical Hall effect. Note also that 
the simple uniform solution above suggests robustness 
against disorder, another element of the integer quantum 
Hall effect. As the transverse gradient of V is indepen- 
dent of g in the uniform bulk, suppose that g is allowed 
to vary slowly from vertex to vertex, more slowly than 
the current. In this regime we would expect quantization 
to persist, and indeed we will confirm that numerically 
later in this paper. 

We will substantiate the assumptions and findings 
above numerically in the next section. 



IV. NUMERICAL RESULTS 

In this section we shall present numerical results for the 
voltages and currents in the wires of the network. These 
numerical studies serve first as a check of the analyti- 
cal results presented in the previous section |III| for the 
case where all the interaction parameters are the same 
for all wires. Second, and more importantly, they serve 
as a stepping stone to the case of non-homogeneous (dis- 
ordered) interaction parameters in the wires, which will 
be considered in Sec. El The method used to solve for 
the voltages and currents in the grid is presented in Ap- 
pendix [Aj 

Let us focus on the armchair layout of Fig. [Ik (simi- 
lar results follow in the case of the zigzag case). Also, 
without loss of generality, we consider below only the x+ 
fixed point. Current is injected and collected uniformly 
into the wires on the right and on the left of the network, 
respectively. More precisely, there are r wires serving 
as connections to the outside on each side of the grid, 
and current / = Ix/r is injected in and collected out of 
these external wires. The total current flowing along the 



horizontal or x-direction is therefore Ix- 

The distribution of the currents in the inner parts of 
the grid that follow from this uniform injection of exter- 
nal currents is shown in Fig. [3j We find a close to uni- 
form distribution, with slightly larger currents closer to 
the edges. This distribution is independent of the value 
of g. These patterns of current flow in the inner wires of 
the grid are in agreement with the current distributions 
discussed in the analytical studies of the previous section. 
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FIG. 3: Currents flowing through the Y-junctions that lie 
along a vertical line in the middle of the bar (x = Lx/2) as a, 
function of vertical position y/Ly. Note that for y values away 
from the edges the currents tend to /i = 1 and 12,3 = 1/2, as 
predicted analytically for the asymptotic limit. 



The Hall voltage is the potential drop Vy along the ver- 
tical or ^/-direction. We note that the potential drop Vy 
is computed by looking at the potentials for two points at 
the same horizontal position (i.e., the same x position), 
one at the top and one at the bottom of the network. 

We show in Fig. |4]the potentials measured at the top 
and at bottom of the (rectangular shaped) grid. No- 
tice that the potentials drop linearly with the horizontal 
direction, but that the difference between the two poten- 
tials, Vy^ is constant. 

The Hall resistance is computed as follows. Let Vy be 
the average over the horizontal positions x of the Hall 
voltage drop. (Since in this case without disorder Vy 
is constant, the average is actually unnecessary here.) 
Then the Hall resistance is given by Rxy = Vy/Ix- We 
find numerically that Rxy = 1, as expected from the 
analytical arguments. Recall that we are working in units 
where e^//i = 1, so indeed we have 



Rx 



h 



(4.1) 



which we find is independent of the value of g. We remark 
that we find that this quantization holds independent of 
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the aspect ratio, orientation (armchair vs. zigzag) or size 
of the grid. 

We also computed the potential difference between 
points on the left and on the right sides of the grid, Vx^ 
as a function of the vertical direction y. In this case 
we find that the horizontal potential difference is almost 
constant as function of y (as opposed to the case of the 
vertical drop Vy^ which is exactly independent of x). The 
difference is bigger, by an amount of order l/Ly^ when y 
is in the middle of the grid as compared to when y is at 
the edges. We define Vx as the ^/-position averaged volt- 
age difference between the left and right sides of the grid. 
The longitudinal resistance is given by Rxx = Vx/Ix^ and 
the longitudinal resistivity by pxx = Ly/Lx Vx/Ix- 

We find that the longitudinal resistance is non-zero, in 



3.5 



agreement with Sec. Ill We find numerically, however, 
that there are finite system size corrections to the ana- 
lytical predictions. We find that 



V3 



g Ly — A(Lx, Ly) 



(4.2) 



where A is a factor of order 1 that corrects for finite sizes. 
We find numerically that in the thermodynamic limit 
^4 ^ 1 for the armchair configuration, whereas A = 
independent of system size in the zigzag case. Therefore, 
in the thermodynamic limit we obtain 



Px 



lim 



Lx 



V3 

9 



(4.3) 



in agreement with the result in Sec. Ill 

The Hall angle Oh is given by idiiiOn = pxy/ Pxx-, and 
we naturally find, given the agreement with the results 
for pxx and pxy above, that 



(4.4) 



in the thermodynamic limit. This Hall angle can be vi- 
sualized very naturally by plotting the voltages at the 
wires on the grid, as shown in Fig. |5] The Hall angle ap- 
pears as the slope of the lines of constant voltage. These 
equipotential lines are straight in this example where all 
the wires have the same interaction parameter g] this is 
no longer the case when disorder is introduced in the next 
SecEl 



bottom edge 




FIG. 4: Voltages at the top and and bottom edges as function 
of horizontal position x/Lx when the node at the top right 
corner is grounded. The grid size is r = 50, c = 60 and 
g — v^- Notice that the difference between the voltages at 
the top and bottom edges for a given x/Lx is exactly 1 in 
natural units. 




FIG. 5: Density plot for the voltages on the grid nodes, for a 
system with r = 50, c = 60 and g — a/3- Notice the constant 
slope of the equipotential lines, which is related to the Hall 
angle Oh • The Hall angle depends on the interaction strength 
and is given by Eq. |4.4[ 



ROBUSTNESS AGAINST DISORDER 



In this section we will generalize the wire networks to 
the case when the interaction parameters g for each of 
the wires in the network are not uniform, but instead 
are drawn independently from a distribution. We shall 
consider a distribution in which g in each of the wires in 
the network takes a value between {g — Sg^g -\- Sg), with 
uniform probability. Because the interaction parameter 
should be positive, 5g < g. 



When the wires connecting to a given Y-j unction have 
different values of ^, the conductance tensor Gjj for a 
chiral fixed point is no longer given by Eq. (2.4), but 
instead it takes the form (see Ref. 25) 



919293 -^91^92+93 



(5.1) 

Using this conductance tensor, one can compute numeri- 
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cally (using the method of Appendix [A| the voltages and 
currents in all wires of the network for a given realization 
of the disorder. 

We shall show below that the quantization R^y = 1 
of the Hall conductance that we found in the clean limit 
remains , in the thermodynamic limit, even in the pres- 
ence of disorder. For a finite lattice, as one should ex- 
pect, there are fiuct nations that we quantify below for 
the armchair configuration. 

We compute Hall resistance Rxy (defined as the aver- 
age of the voltage differences between top and bottom of 
the network, divided by the injected current) for several 
realizations of disorder and system sizes. For a fixed sys- 
tem size, we then find the disorder average Rxy and stan- 
dard deviation AR^y = \J R^y — Rxy^ of Rxy We find 
that Rxy ^ 1 as the number of realizations increase, and 
that the standard deviation ARxy ^ as L increases (we 
use lattices with r = c = L). We show in Fig. |6]the finite 
size scaling of the ARxy. That ARxy ^ in the thermo- 
dynamic limit means that the system is self- averaging, 
and therefore Rxy 1 independent of disorder in the 
thermodynamic limit. We conclude then that quantiza- 
tion is robust against disorder. 
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FIG. 6: Standard deviation of the Hall resistance for 100 
simulations with g — and = ^/lO as a function of 1/L 
for a grid with r — c — L (which fixes the aspect ratio). It 
scales to zero in the large L limit, implying that the system is 
self- averaging and the Hall resistance Rxy 1 independent 
of disorder in the thermodynamic limit. 

We have also checked the effects of disorder for the 
zigzag configuration, reaching similar conclusions that 
disorder does not alter the quantization of the conduc- 
tance in the thermodynamic limit. 

In summary, we find that, in the thermodynamic limit, 
the general results of the previous sections hold even in 
the presence of disorder. 



VI. CONCLUSIONS 

We investigated the transport properties of hexagonal 
networks whose nodes are Y-j unctions of quantum wires. 
In our model the 3x3 conductance tensor for each Y- 
j unction is dictated by the low energy RG fixed point, 
but the transport is treated classically between any two 
junctions. We find a surprising result: in spite of relax- 
ing quantum coherence between the junctions, we find a 
quantized Hall resistance. 

Specifically, in the regime where the junction conduc- 
tance is controlled by the chiral fixed points x± (when 
the interaction parameter obeys 1 < ^ < 3), the network 
exhibits a quantized Hall resistance: Rxy = ±/i/e^. This 
quantization is similar to that in the integer quantum 
Hall effect. Further, the quantization is independent of 
the interaction parameter g even in the presence of disor- 
der in g. The quantization of the Hall resistance follows 
from the specific form of the conductance tensor at the 
RG stable chiral fixed point at each Y-junction. However, 
unlike in the quantized Hall effect, where the longitudinal 
resistivity vanishes, pxx is not zero: pxx = {V^/9)h/e^. 
Dissipation in the longitudinal direction is a result of de- 
coherence within the wires. We emphasize that in our 
model the wires are classical, but the nodes remain quan- 
tum mechanical and the form of the conductance tensor 
G at each junction is constrained by quantum scattering 
effects. 

Finally, let us comment on the finite temperature cor- 
rections to the value Rxy = ±/i/e^ in the network. As 
opposed to the case of the quantum Hall effect where the 
quantization is exponentially accurate because of an en- 
ergy gap, the quantization in the networks has a power 
law correction in T because the wire networks are gapless. 
The quantization should be as accurate as the conduc- 
tance tensor is close to that of the RG fixed point. The 
corrections to the conductance tensor scale as T^, where 
A = 4^/(3 + ^^) is the scaling dimension of the leading 
irrelevant operator at the chiral fixed points.^ ^ 
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Appendix A: Method 

Our numerical approach consists of solving the full 
lattice model exactly. In this section we describe our 
methodology in detail. 

Consider an arbitrary lattice with r external wires on 
each side and 2c external wires at the top and bottom 
edges. An equal current will be injected into each of the 
r wires on the left, and the 2c edge wires will have a 
current of zero. This defines the boundary conditions. A 
2x2 lattice is shown for example in Figure jT] For later 
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convenience we include a row of 2c "ghost nodes" , shown 
as dotted hnes at the top edge, but they are only there 
to facilitate the numbering scheme and no current will 
flow through them. Including the ghost nodes there are 
a total of N = Ac{r + 1) nodes. 



(I2>i — Hl3 



(loh — Hii) 



(2d}j — H21 



(I8h — Hl9) 




FIG. 7: Example with r = 2, c = 2. Note the row of "ghost 
nodes" at the top edge. 



The points on each wire that emanate from each node 
are governed by the equation V = GI where G is a 
3x3 matrix. Thus we start with 3A/' degrees of freedom. 
However, starting in this way introduces many redun- 
dant variables in the bulk because in a classical wire the 
current is the same everywhere along the wire and so is 
the potential. We will unify the two points on each wire 
in the bulk by imposing a set of constraints. In general 
there are 6cr + c — r such constraints, which equals the 
number of wires in the bulk. 

To write down the full network equation let us label 
each of the 3A/' points by (n, z), where n = 1, . . . , A/" is 
the node index and i = 1,2,3 refers to the point on each 
wire that emanates from each node. The potentials and 
currents at each of these points are denoted by V^"^^ and 

, respectively. To illustrate this notation, in Fig. [t] 
the constraint along the wire that connects nodes 6 and 
7 would be written as V^^^ = v}^^ and /f ^ = -/f \ 

Each node obeys the relation y(^ ) — qMjM where 
G^^^ is the 3 X 3 matrix from Eq. The network 



/(2) G(2) 



\ (v^^^\ 



cm J 



(Al) 



Next we impose the constraints to reduce the effec- 
tive dimensionality of the problem. Start with the set of 
point pairs on each wire in the bulk {(n, i), (m, i)}, where 
n and m are nearest neighbor nodes. The constraints are 
f/^^ = y/^^ and = — /|^^ for each pair. We im- 
pose the constraint on voltages by adding the (n, i)-th. 
and (m, i)-th columns together, removing the (m, i)-th 
column and removing 1^^^^ from the vector of potentials 
in Eq. ( |A1| ). Similarly we impose the constraint on cur- 
rents by adding the (n, i)-th and (m, i)-th rows^ deleting 
the (m, i)-th. row and removing from the vector of 
currents. Also we replace the current I^^^ that has not 

been eliminated by zero because 4''^ + 4"^^ = 0. There- 
fore each constraint is equivalent to removing one row 
and one column and reduces the dimensionality of the 
original problem by one. Furthermore, we have replaced 
each current in the bulk by zero which is important be- 
cause the only currents that are left in Eq. ( Al ) are fully 



determined, being equal to either zero in the bulk or to 
the boundary conditions. 

Eliminating the ghost nodes is straightforward - we 
simply remove the ghost currents, potentials and their 
associated rows and columns in Eq. (Al). This reduces 



is thus described by the following linear equation with 
constraints: 



the dimensionality further by 3 x 2c, which is the number 
of wires emanating from the ghost nodes. The flnal step 
is to flx the gauge. Since all potentials are determined 
up to an overall constant, we pick an arbitrary potential, 
set it to z ero, and remove the associated row and column 
from Eq. ( Al ). 

To summarize, we started with 37V = 12c(r + 1) re- 
dundant degrees of freedom and then through successive 
transformations we imposed 6cr + c — r constraints in the 
bulk, eliminated 6c ghost points, and flxed one poten- 
tial to zero. The dimensionality has thus been reduced 
to 6rc + 5c + r — 1 and, crucially, the only currents ap- 
pearing are either zero or flxed by boundary conditions. 
Having eliminated all redundancies allows us to solve for 
the potential at any point, as a function of the boundary 
currents, by inverting the reduced version of Eq. (Al), 
which we do numerically. 

The generalization to random couplings g is straight- 
forward. The derivation proceeds in exactly the same 
way as we just described, but we start with non- uniform 
G(^). 
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